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DGDcription 

METHOD FOR IDENTIFICATION OF AN OSCILLATION IN AN 
ELECTRICAL POWER SUPPLY SYSTEM 

CLAIM FOR PRIORITY 

This application claims priority to International 

Application No. PCT/DEOO/03434 which was filed in the 
German language on September 21, 2 000. 

TECHNICAL FIELD OF THE INVENTION 
The invention relates to a method for producing at 
least one signal (oscillation signal) , which indicates 
an oscillation in an electrical power supply system-; — ±^ 
which — method — fcfee — phaoc — current — and — t-he — phase — voltage 
arc in each caoc oamplcd from at — Icaot one phaoc of — t-he 

power — supply — oystcm, forming — phaoc — current — aftd — phaoc 

voltage oamplc valuco, — impcdancQ valuco — are formed from 
the phaoc — current — and phaoc voltage — oamplc valuco^ — feiie 
impedance valuco — a^ee — monitored — §ei= — fcfee — prcoGnce — — aay 
oocillation — and, — — etn — oocillation — is — identified, — afe- 
Icaot — efie — memory — clement — is — oot , — afid — fefee — oocillation 

oignal — is — output — &^ — ifes — output , after — oetting — fefee 

memory clement, — further impedance valuco arc Gheckcd to 
dot ermine — whether — fefee — oocillation that hao — been found 

4r& otill continuing, fe4=ie memory element rcmaino 

uninfluenced — i^ fehe oocillation — continueo, aftd — 

memory clement io reoct if the oocillation hao otopped . 



A method ouch ao — thio — io dcocribcd in German Laid- 
Open Specification DE 195 03 626 Al describes a method 
of identifying an oscillation . In this method, once a 
memory element has been set, further impedance values 
are checked to determine whether the oscillation that 
has been found is still continuing, by determining the 
rate of change of the magnitude of respectively 
successive impedance values — and, — i #. If it is found 
that the rate of change is above a limit value, this 
identifies that the oscillation has stopped, and the 



BACKGROUND OF THE INVENTION 



memory element is reset. It has — bocn — found — t-e — fee — 
difficult to define such a limit value, particularly 
when a large number of generators are connected in the 
power supply networks, and complex oscillations can 
5 thus occur. 

SUMMARY OF THE INVENTION 
The invention relates to a method for producing at 
least one signal (oscillation signal) , which indicates 

10 an oscillation in an electrical power supply system, in 
which method the phase current and the phase voltage 
are in each case sampled from at least one phase of the 
power supply system, forming phase current and phase 
voltage sample values, impedance values are formed from 

15 the phase current and phase voltage sample values, the 
impedance values are monitored for the presence of any 
oscillation and, if an oscillation is identified, at 
least one memory element is set, and the oscillation 
signal is output at its output, after setting the 

2 0 memory element, further impedance values are checked to 

determine whether the oscillation that has been found 
is still continuing, the memory element remains 
uninfluenced if the oscillation continues, and the 

- - memory el-ement is reset if the- oscillation-has stopped, 

25 

The invent ion is — baood — or — fehe — obj act — ef^ — opacifying 
specifies a method of detecting fey — mcana — ef — which — the 
oscillation behavior of an electrical power supply 
system can be dctGctcd oafcly and reliably at — aii — timco 

3 0 in a safe and reliable manner . 

Fe^^ — a — method — ef — the — type — mentioned — initially, — thio 
obj GCt — i^^ — achieved — according — fe-e — In one embodiment of 
the invention , there is a method in that tho checks e# 
3 5 the further — impedance values and makes use of an 
oscillation model which is formed from previous 
impedance values associated with the oscillation, or 
from variables which are dependent on these impedance 
valuco ; avalues . A check is then carried out to 



determine whether a further impedance value formed at 
that time or a variable which is dependent on thio the 
further impedance value differs from the oscillation 
model, and any occurrence of a further impedance value 
5 or of a variable dependent on this impedance value 
which differs from the oscillation model is assessed as 
the oscillation having stopped, 

One major One advantage of the method according to the 
10 invention is that the oscillation model allows even 
complex oscillations to be described, and it is thus 
possible to identify that the oscillation has stopped 
with a high level of reliability even in the case of 
such complex oscillations. 

15 

The oscillation model can advantageously be 

determined by means of a least squares estimation 
method. This estimation method allows a mathematical 
oscillation model to be produced from successive 
20 impedance values which have been formed after the 
setting of the memory element, that is to say after the 
start of the oscillation. 

' A "function in the form f (x) = ax^+bx^+cx+d with the 
25 parameters a, b, c and d can be used as the model rule 
for thio the oscillation model, in which one or more 
parameters can be defined to be zero before the start 
of the estimation method. First, second or third order 
power functions ca nthuo be used as the model rule. 
30 Furthermore, a sum of sine and cosine functions, which 
decay with time, can be used as the model rule for the 
oscillation model- These model rules make it possible 
to describe even complex oscillations mathematically. 

35 The oscillation model can be formed directly for 

the determined impedance values of the oscillation, or 
else for variables dependent on these impedance values. 
Resistance values R, reactance values X, time 
derivative values dZ/dt of the impedance, time 



derivative values dR/dt of a resistance or time 
derivative values dX/dt of a reactance can be used as 
dependent variables. Choice of the most suitable 
variable for the oscillation model makes it possible to 
5 determine with a high level of reliability that the 
oscillation has stopped, with the choice of the 
variable being dependent on the individual system 
configuration of the electrical power supply system. 

10 In one advantageous embodiment of the invention, 

positive phase sequence system impedance values can be 
formed from the phase current and phase voltage sample 
values, and a common memory element can be provided, 
and a common oscillation signal can be produced, f or all 

15 the phases in the power supply system. This variant can 
be used when the aim is to make a statement relating to 
any oscillation occurring at the same time in all the 
phases in the power supply system. 

2 0 In a f urthe radvQntagcouo embodiment of the method 

according to the invention, phase impedance values are 
formed from the phase current and phase voltage sample 
values of each phase of the power supply system to be 
investigated for oscillation, and- a dedicated memory 

25 element is provided, and a dedicated oscillation signal 
is produced, for each of these phases. In this 
embodiment, the oscillation response of each individual 
phase in the power supply system can be investigated 
ocparatoly, that — is — fee — say separately . That is, both 

30 the starting and the stopping of an oscillation are 
identified. This is particularly advantageous when 
oscillations occur i nonly a single phase, but not in 
all the phases, in the power supply system. 
Oscillations such as these frequently occur in the case 

35 of so-called single-pole pauses in high-voltage 
systems. Single-pole pauses are produced by single-pole 
conductor ground faults, which can be expected 
frequently in high-voltage systems, and in which an arc 
is struck between one conductor and ground. In this 



type of fault, a single-pole pause is produaod — that — ie 
^ -Ba yproduced . That i s , the single phase in which the 
single-pole conductor -ground fault has occurred is 
switched off briefly. The arc is thus quenched, and the 
5 fault is frequently corrected. Switching off a single 
pole of one phase can result in oscillations occurring 
in the remaining phases which are not switched off. 
These oscillations cannot be identified, for example, 
by monitoring the positive phase sequence system 

10 impedance values, since positive phase sequence system 
impedance values can — b eare f ormed only when sample 
values are available f or all the phases in the power 
supply system. In the case of a single -pole pause, it 
isftew — highly advantageous to be able to produce a 

15 dedicated oscillation signal for each phase in the 
power supply oyotom; — dur ing system , During the single - 
pole pause, this oscillation signal is producedeaiy for 
those phases which are not switched off. The 
oscillation behavior of the power supply system can 

2 0 thus be determined individually for each phase, and 

independently of the state of the other phases. 

The phase impedance values of the individual 

- - phases in the electrical power supply system can, -for 
25 example, be formed by, in order to form the phase 

impedance values , 

a variable U_re containing including the real part 

of the phase voltage sample values, a variable U_im 

containing including the imaginary part of the phase 

3 0 voltage sample values, a variable 

containing including the real part of the phase current 
sample values and a variable I_im containing including 
the imaginary part of the phase current sample values 
being formed from the phase current and phase voltage 
35 sample values (i, u) for each phase, 

a phase real power variable P being determined 
from P = U_re • I_re - U_im"I_im^ 

a phase wattless component variable Q being 
determined from Q = U_im-I_re + U_im-I_re, 



a squared phase current amplitude variable 
being determined from = I_re • I_re + I_im-I_im^ 
system- frequency components in each case being 
removed by means of least squares estimation 
5 method from the phase real power variable P, from 

the phase wattless component variable Q and from 
the squared phase current amplitude variable 1^ , 
and 

phase resistance values R being determined from 
10 R=P/I^ and phase reactance values X being 

determined from X=Q/1^, and phase impedance values 
Z=R+jX being thus determined. 

When forming the phase impedance values, it 

15 is part icularly advantageous to remove the system 
frequency components (for example 50 Hz components) 
from the phase real power variable P, from the phase 
wattless component variable Q and from the squared 
phase current amplitude variable by means a least 

20 squares estimation method in each case. Such system 
frequency interference components would adversely 
affect the evaluation of the phase impedance values 
determined from these variables . 

25 BRIEF DESCRIPTION OF THE DRAWINGS 

In order to explain the invention further. 

Figure 1 shows a block diagram of an exemplary 
embodiment of the method according to the invention. 
invention, — Figure 2 Figure 2 shows a block diagram 

30 for determining the phase impedance valuco, — values . 

Figure 3 shows the impulse response of the filters 

used for impedance determination, determination. 

35 Figure 4 shows the real power and wattless 

component variable profiles before filtering, 
aft ^f iltering . 



Figure 5 shows the real power and wattless 




component variable profiles after filtering. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

Figure 1 shows, schematically, a method for 

5 determining the oscillation behavior of a three-phase 
electrical power supply system, by means of which a 
dedicated oscillation signal Pdl, Pd2 and Pd3 is 
produced for each phase in the power supply system. 
This is done by providing three changeover switches Ul, 

10 U2 and U3 and three memory elements Spi, Sp2 and Sp3 ; 
thc Sp3 , The connecting lines between the individual 
units in the layout are designed with three poles. The 
phase current and phase voltage sample values i and u 
of all three phases are supplied to a unit for 

15 impedance determination lb, at whose output phase 
impedance values Z are output for the three phases . 
These phase impedance values Z are supplied via the 
changeover switches Ul, U2 and U3 to an oscillation 
identification unit Pe . The oscillation identification 

2 0 unit Pe uses the time profiles of the phase impedance 

values Z to identify the occurrence of any oscillation 
in the individual phases, for example in phase 1, and 
emits an oscillation set signal Ps at its output for 
each phase in which oscillation is identified, for 
25 example for the phase 1. The oscillation set signal Ps 
sets the memory element associated with the respective 
phase, for example Spl, and this memory element emits 
at its output the phase-specific oscillation signal, 
for example, Pdl . When an oscillation signal is being 

3 0 emitted, for example the oscillation signal Pdl, the 

changeover switch, for example Ul, associated with the 
respective phase is switched over. The phase impedance 
values Z which are still formed for the phase in which 
the oscillation has been identified, for example the 
35 phase 1, are supplied to an oscillation signal 
resetting unit Pii . This oscillation signal resetting 
unit Pvi identifies that the oscillation has stopped 
and, in this case, emits an oscillation reset signal Pr 
at its output, which resets the memory element for the 



respective phase, for example Spl . The oscillation 
signal for the respective phase, for example Pdl, is 
thus also no longer emitted, and the respective 
changeover switch, for example Ul, moves back to its 
5 original position once again. In response to a 
stimulus, a unit for phase selection Pa ensures that 
the phase impedance values of the phases to be 
investigated for oscillation are in each case processed 
by the oscillation identification unit Pe and by the 
10 oscillation signal resetting unit Pii. 

The method of operation of the four units 

comprising impedance determination lb, the oscillation 
identification unit Pe, the oscillation signal 
15 resetting unit Pii and phase selection Pa will be 
explained in more detail in the following text. 

^As shown in Figure 2 , the phase current and phase 

voltage sample values i and u are filtered in the 

20 impedance determination unit lb by means of orthogonal 
FIR filters Fl, F2 , F3 and F4,. thus resulting in the 
production of a variable U_re containing the real part 
of the phase voltage sample values, a variable U__im 
containing the imaginary part of the - phase- voltage 

25 sample values, a variable I_re containing the real part 
of the phase current sample values, and the variable 
I_im containing the imaginary part of the phase current 
sample values . 

30 Figure 3 shows the impulse responses of the 

filters Fl to F4, with the impulse response of the 
filters Fl and F3 which determine the real parts being 
annotated "o", and the impulse response of the filters 
F2 and F4 which determine the imaginary parts, being 

35 annotated "+". 

^As shown in Figure 2, following this, a phase real 

power variable P is calculated in accordance with 
equation (1) below and a phase wattless component 



# m 

variable Q is calculated in accordance with equation 

(2) in the unit 5, and a squared phase current 
amplitude variable is calculated in accordance with 
equation (3) in the unit 6. 

P = U_re • I_re - U_im-I_um (1) 

Q = U_im-I_re + U_im-I_re (2) 

= I re • I re . + I im-I im (3) 



^After this, the phase real power variable P, the 

phase wattless component variable Q and the squared 
phase current amplitude variable are filtered in 

10 units 7 and 8 in order to remove the interference 50 Hz 
components contained in these variables; this results 
in the filtered variables P', Q' and I^'. The least 
squares estimation method used for this filtering will 
be explained in detail further below. 

15 

The upper illustration a) in Figure 4 shows the 

profile of the real power variable P, and the lower 
illustration b) shows the profile of the Wattless 
component variable Q, before filtering by means of the 

20 least squares estimation method, in each case plotted 

_ - .against -the time t. - 

The upper illustration a) in Figure 5 shows the 

profile of the real power variable P', and the lower 
25 illustration b) shows the profile of the Wattless 
component variable Q', after filtering by means of the 

least squares estimation method; jrfe rmethod- Tt can 

clearly be seen that the 50 Hz components have been 
removed . 

30 

^As shown in Figure 2, after the filtering in the 

unit 9, phase resistance values R and phase reactance 
values X are determined in accordance with equation 
(4) , and the phase impedance values Z=R+jX determined 
35 from them are emitted at the output of the impedance 
determination unit lb. 



10 



R=pi/l2i X=Q*/I^' (4) 

least squares estimation method using a signal 

model in accordance with equation (5) is applied 
separately to each of the variables P, Q and in 
order to filter out the 50 Hz components 
containcd included in the phase real power variable P, 
in the phase Wattless component variable Q and in the 
squared phase current amplitude variable I^. 

y^=A-e ^ •sin(cDok T^) + B e"^ •cos(cOok TA) + C (5) 



The estimation method uses the squared phase 

current amplitude variable I^, the phase real power 

15 variable P and the phase Wattless component variable Q 
to calculate the parameters A, B and C in the signal 
model . The parameter C gives the sought magnitude of 
the phase real power variable P', of the phase Wattless 
component variable Q* and of the squared phase current 

20 amplitude variable I^*. The summands with the 
parameters A and B model the 50 Hz components. The 
-variable - (DO - is the frequency (50 Hz) -to be- filtered 
out, and Ta is the sampling time. 

25 If an equivalent circuit with only two generator 

machines at the ends of a power transmission line is 
considered for the power supply system, the amplitude 
of the 50 Hz components decays with the time constants 
T of the sum impedance between the two generator 

3 0 machines in accordance with equation (6) , where L is 
the loop inductance and R is the loop resistance of the 
circuit which is closed via the two generator machines. 



r = 



2£ 



35 



(6) 



The coefficients A, B and C are determined such 

that the sum of the squares of the errors between 
values y determined from the phase current and phase 
voltage sample values i and u and the sample values yk 
calculated in accordance with equation (5) becomes a 
minimum (see equation (7)). 



J= X(yi-h(e,))'-^MiN 



(7) 



10 



15 



In equation (7) , J represents the Q-criterion to 

be minimized. The signal model included in equation (5) 
is used as the function h(©k) . The parameters A, B and 
C to be determined form a vector ©j^ in accordance with 
equation (8) . 



B 



(8) 



The Q-criterion 



IS 



derived based 



on 



the 



parameter vector in order to solve the minimization 
2 0 task. For the signal model in accordance with equation 
(5) , this then results in equation (9) together with 
equation (10) . 



0= S2Y;^(yi-LQ0 



(9) 



cos^ 



(10) 



If equation (9) is solved for the vector 0^/ then 

this results in equation (11) , by means of which, and 
using the matrix Sk included in equations (12) and 
(13) , the vector ©jc is determined. 



. f2n 
sin 



(11) 



(12) 



e ^ 



1 



(13) 



10 



Of the parameters A, B and C GontainGd included in 
only the parameter C is evaluated. The 



the vector Qv 



vectors 



in accordance with equation (10) and the 



matrix S^k in accordance with equation (13) 
calculated and are stored as constants, so that 
are available every timc when the method is used. 



are 
they 



15 ^Monotony criteria are applied to the locus curves 

of the impedance values in the impedance plane in the 
oscillation identification unit Pe, in order to 
identify the oscillation process. This method for 
identification of the oscillation process is known — pe*^ 

20 ee-? — and is described in German Patent DE 197 46 719 CI. 



The oscillation signal resetting unit 



Pii 



determines whether an oscillation which has already 
been identified is still continuing. The procedure used 
25 for this purpose comprises the production of an 
oscillation model for phase impedance values Z 
associated with the oscillation. A check is then 
carried out to determine whether the locus curve which 
is described by the newly determined phase impedance 



10 



♦ 



values Z still corresponds to the oscillation model. 
When producing the oscillation model, it is assumed 
that the locus curve is free of discontinuities, and 
that its direction changes only — very slowly. In the 
present exemplary embodiment, the locus curve is 
described by a first order power function, that is to 
say by a linear equation, in accordance with equation 
(14) . 

X(R)=m-R+Xo (14) 
The parameters m and XO are determined by means of 



a non-recursive least squares estimation method from 
the last N determined phase impedance values Z. 

15 The linear equation is used as a model rule for 

the least squares estimation method, with the parameter 
m characterizing the gradient, and the parameter XO 
characterizing the offset of the linear equation. The 
parameters m and XO for the model in accordance with 

2 0 equation (14) are determined from the last determined 
value pairs (Ri, Xi) of the phase impedance values Zi 
such that the sum of the squares of the errors between 

„ . the . _ values Xi ^ determined„ _f rom ,the measured _ phase . 

current and phase voltage sample values i and u and the 

25 values X calculated in accordance with equation (14) is 
a minimum (see equation (15) ) . 



1 



E(^.^h(e,F->MIN (15) 



30 In equation (15) , J is the Q-criterion to be 

minimized, the model rule in accordance with equation 
(14) is used as the function h(0k). In accordance with 
equation (16) , the parameter vector 0jc contains the 
parameters m and XO, to be determined, from the model 

3 5 rule. 




(16) 



In order to solve the minimization task, the Q- 

criterion J muat — feeis derived on the basis of the 
5 parameter vector ©k- The following equations (17) and 
(18) are then obtained for the signal model in 
accordance with equation (14) : 



10 



0= t^r]{x,"h{e,X) (17) 



(18) 



If equation (17) is solved for the parameter 

vector ©k, then this results in equation (19) for 

determining the parameter vector ©k. 



15 



where 



(19) 



(20) 



(21) 



^After substitution of the parameters in equation 

20 (14) , this thuo results in the estimated oscillation 
model. If a newly determined phase impedance value Z 

corresponds to the oscillation model , that i-s fee 

e-a ymodel . That is^ it is within a tolerance band 

around the linear equation represented by equation 
25 (14) , then this identifies that the oscillation is 
continuing. If the newly determined phase impedance 
value Z is outside the tolerance band, then this 
identifies that the oscillation has stopped, and an 
oscillation reset signal Pr is emitted at the output of 



4 4 

the oscillation signal resetting unit Pii for the memory 
element Spl, Sp2 or Sp3 for the respective phase. 

The phase selection unit Pa receives a stimulus 

5 from, for example, a distance protection device which 
is not illustrated. Depending on the nature of the 
stimulated loop, it determines the phases for which the 
oscillation identification unit Pe and/or the 
oscillation signal resetting unit Pii should investigate 
10 the oscillation behavior. The following table shows the 
association between the stimulated loops and the 
phases . 



Stimulated loops 


Phases to be investigated 
for oscillation behavior 


LIE 


LI 


L2E 


L2 


L3E 


L3 


L12 


LI and L2 


L23 


L2 and L3 


L31 


LI and L3 



